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Abstract 



> 

We study the transition of a heavy quark pair from octet to singlet color configurations at 
next-to-next-to-leading order (NNLO) in heavy quarkonium production. We show that the 
infrared singularities in this process are consistent with NRQCD factorization to all orders 
in the heavy quark relative velocity v. This factorization requires the gauge-completed 
matrix elements that we introduced previously to prove NNLO factorization to order v 2 . 



a: 

1 Introduction 



Heavy quarkonium production serves as a testing ground for perturbative and effective field 
theory treatments of QCD, particularly nonrelativistic QCD (NRQCD) [T]. NRQCD, which 
relies on an expansion in the pair relative velocity, as well as in a s , has provided compelling 
explanations for quarkonium production at collider [2j and fixed target experiments At the 
same time, puzzles remain, especially from polarization measurements at the Tevatron [3], and 
associated production at the B factories [H]. Heavy quarkonium production in the evolving QCD 
medium at the RHIC and LHC jB] may also play a role in the detection and analysis of new 
states of strongly-interacting matter. For recent updates on heavy quarkonium production at 
zero and finite temperature QCD see Refs. IE] • 

The application of NRQCD to heavy quarkonium production processes is based on a very 
specific factorization property [I], for which a complete proof has not yet been developed [B]. With 
this in mind, we tested the factorization hypothesis at next-to-next-to-leading order (NNLO) in 
Ref. We found that, as a necessary condition for factorization, the conventional 0(v 2 ) octet 
NRQCD production matrix elements must be redefined by incorporating Wilson lines that make 
them manifestly gauge invariant. We referred to these as gauge-completed matrix elements. 
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This result was derived by employing an eikonal approximation for the coupling of soft radi- 
ation to the heavy quarks. To order v 2 , this approximation is adequate to treat the lowest-order 
electric dipole transitions ^Oj that transform an octet pair to a singlet, without modifying the 
spin. In this paper, we extend our eikonal NNLO analysis, and show that the factorization of such 
transitions can be extended from order v 2 to finite v, that is, to all orders in the relative velocity, 
in terms of gauge-completed NRQCD matrix elements. The essential result is that at all orders 
in v, the infrared divergence at NNLO is independent of the direction of the light-like Wilson 
line that renders the matrix element gauge invariant. We find this result intriguing, and while 
this modest extension of our previous result does not address the behavior of spin-dependent 
operators, it should encourage further work on the factorization theorem. 

In the following section we review the role of gauge-completed matrix elements in NRQCD 
factorization, and the requirements that factorization places on these matrix elements and the 
infrared poles in dimensional regularization that they must match. We also introduce the specific 
NNLO eikonal factor that we will calculate to test factorization at this order, and give the result 
of our calculations. The details of our NNLO calculation are presented in Sec. 3, where we verify 
that the necessary conditions for NRQCD factorization are met to all powers in the relative 
velocity, after which we give a brief conclusion. 

2 Infrared Poles and NRQCD Factorization 

In NRQCD, the production cross section for heavy quarkonium H at transverse momentum px 
factorizes into a sum of perturbative functions times universal matrix elements, 

da A +B^H+x(PT) = ^2 d °A+B^cc[n]+x{PT) {O") , (1) 
n 

where each matrix element (O^) represents the probability for a heavy quark pair in state [n] to 
produce quarkonium state H. The states n may, in particular, be color octet or singlet. 

Correspondingly, at large Pt, the fragmentation function for parton i to evolve into a heavy 
quarkonium is factorized according to [TU] 

D H /i(z,m c , /i) = 2^ di->cc[n\{z, m c , n) (O") , (2) 

n 

in terms of the same matrix elements, along with perturbative functions di^ c5 [ n ](z,m c , fi) that 
describe the evolution of an off-shell parton into a quark pair in state [n], including logarithms 
of \ijm c . 

Although we cannot compute the full fragmentation function in perturbation theory, for 
NRQCD factorization to hold we must be able to compute the function d^^Az, m c , fi) in 
Eq. (J2J) systematically. To do this, we must be able to use regularized perturbation theory to 
compute probabilities for the production of a singlet quark pair from any local source, and match 
the infrared singularities of these probabilities to the matrix elements in the sum over n in Eq. ([TJ 
and/or Eq. (0). 
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2.1 Gauge-completed matrix elements 

As advocated in jH], it is natural to define gauge-invariant octet NRQCD operators Off , of the 
general form 

(off(o)) = (o x ] (ok, b v(o) ^! A)t (o) e6 foto H ) $; A) (o) 6a ^(o)< a x(o) \ (3) 



in terms of heavy quark (ip) and antiquark (x) operators, and local combinations of color and spin 
matrices and/or covariant derivatives, denoted by n n ^ e and K' na . The operator a H creates quarko- 

nium H, and the operators $^(0) are Wilson lines, that is, ordered exponentials, constructed 
from the gauge field in adjoint matrix representation, as 



$J A) (0) = Pexp 



ig / d\l- A (A \l\) 



(4) 



where V denotes path ordering and is the velocity of the source. In (complex conjugate) 
amplitudes, (anti)time-ordering is understood. 

For NRQCD factorization to hold, a necessary, and superficially paradoxical, property of the 
gauge-completed matrix elements is that their long-distance behavior must be independent of 
the vector that we choose to define them Such a dependence would be inconsistent with 
NRQCD factorization, because the infrared divergences of Off must match those of cross sections, 
in which there is no information on P. In Ref. jH], we have verified the /-independence of the 
infrared pole to order v 2 in the relative velocity of the pair, at order NNLO. We will extend this 
result below, to all powers in v, again at NNLO. 

2.2 Matrix elements and infrared universality 

In this paper, we will study the infrared behavior of the octet, S-wave matrix elements 

M^ I \P 1 , P 2 , 0=£ <0| X f (0) l?ty(0) $J A)t (0U I [c{P,)c{P2)] {I) x) 
x 

x(X [c(Pi)c(P 2 )]« | $J A) (0) fea ^(0) T^ X (0) |0> , (5) 



where I = 1, 8 labels the color of the heavy quark pair, with momenta Pi and P 2 . The Tf q ^ are 
color generators in the quark fundamental representation. We will compute the infrared poles of 
Eq. (j2J) to NNLO. At lowest order of course, only 1 = 8 contributes in the final state. At higher 
orders, however, the octet state mixes with the singlet state, through radiation to states X. The 
phase space of X is cut off at a UV scale fi, which we take to be of the order of the heavy quark 
mass. 

In Eq. we specify the pair's momenta in conventional fashion, in terms of the total pair 
(P) and relative momentum (q), and center-of-mass velocity (v) by 

^1 = ^ + 9, p 2 = ^-g, ^ 2 = |^' (6) 
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where 2E* is the total center of mass energy of the pair in the pair's rest frame. We note that 
our convention for the relative velocity here differs from the one in our previous study, Ref. [5], 
where we defined the relative velocity as 2q/m, with m the heavy quark mass. 

NRQCD requires that when we expand M. of Eq. (j^J in powers of q in perturbation theory, 
we should find infrared-finite coefficient functions times infrared- sensitive but universal NRQCD 
matrix elements, 

M^iP^P^l) =J2^ 8) (PuP2J) (Oi) , (7) 

n 

where Cn^ is a perturbative coefficient function for NRQCD operator G^ with the sum over all 
operators n for a given final state in which the heavy quark pair has color /. To test these ideas, 
we must study the infrared behavior of such a matrix element when the color of the final state 
is fixed as a singlet, 1=1. 

The formation of a heavy quarkonium state, of course, cannot be realized to any fixed order 
of perturbation theory. Nevertheless, the factorization expressed in Eqs. (JTJ), (J2J) and (JJJ) is useful 
to the extent that we can systematically compute corrections to short- distance functions in each 
case. This, in turn, requires that the infrared poles encountered in the evolution of the heavy 
quark pair from octet to singlet in (G„) in Eq. (J2J) match those of cross sections for the production 
of a singlet quark pair. As already noted, a necessary condition for this matching is that the poles 
in the matrix element should not depend on the direction of the vector Z M that is introduced in 
gauge completion of the matrix elements, Eq. Q, since the choice of Z M is matter of convention. 

Now we are ready to describe the NNLO calculation that tests these ideas. A full NNLO 
calculation of Ai^ 1 ^ in Eq. (J5J) would be impractical, but its infrared singularities are easier to 
compute. These divergences can be generated by a factorization that is much simpler than Eq. 
(J7J), and which can be carried out at fixed momenta Pi and P 2 , 

M^iPt, P 2 , 1) = J^M^i, P 2 , 1) £^\P U P 2 , s) . (8) 

j 

All spin information in Ai^ 1 ^ is contained in another short- distance function C 8 j, which de- 
scribes a transition of octet to color configuration J at short distances. Although C 8 j may depend 
on Z M , it must be finite for e — ► 0. All l/e poles are absorbed into an infrared factor, S^ 1 ^, 
whose pole structure, however, must be independent of l^. The expansion of the infrared factor, 
£( J ~^ 1 )(P 1 ^ P 2) e) in the relative velocity, v of the pair should lead us back to a set of /-independent 
operator matrix elements for producing a color singlet quark pair (/ = 1) in Eq. (J2J). The essen- 
tial result of this paper is that at NNLO the infrared poles of the function £( 8_>1 '(.Pi, P 2 , s) are 
indeed independent of the vector /, for arbitrary v. This generalizes the result of Ref. [5] from v 2 
"electric dipole" transitions to arbitrary powers of v at NNLO. 

Specifically, we will find an explicit single-pole contribution, which can be written as a pref- 
actor determined by N c , the number of colors, times a function that depends only on the relative 
velocity, 

£W(Pi,P a ,e) = -^(N 2 -l)X^(v,e). (9) 
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The velocity-dependent factor X 8 ^ 1 (f ,e) is given by 



2/(1*1) 



In 



1 + 

i-/(N) 



(10) 



where (u = \v\) is the relative velocity of the quark and antiquark in the pair center of mass, and 
where /(f) = 2v/(l + v 2 ). As anticipated, Z 8- * 1 ^, e) is independent of /. We will compute the 
color factor given in Eq. Q below. To derive the result of Eq. (jlUJI . we will first recall the use of 
the eikonal approximation to isolate infrared behavior. 



2.3 Soft gluon interactions 

The eikonal apporoximation reproduces all infrared divergences in the evolution of the pair into 
the final state, and it must be defined by an infrared regularization. We will use a continuation 
to 4 — 2e dimensions, with e < 0. 

The eikonal approximation for the interactions of the heavy quarks with soft gluons is gener- 
ated by ordered exponentials, this time in fundamental representations and in the directions of 
the heavy quark and antiquark momenta. The perturbation theory rules for the ordered expo- 
nentials are equivalent to the eikonal approximation. Eikonal quark propagators and gluon-quark 
vertices are specified respectively by 



(P-k + ie) 



(11) 



where the plus is for the antiquark and the minus for the quark vertices, and where (3^ is a 
time-like four- velocity. Because the product of an eikonal propagator and vertex is always scale 
invariant, we will use below the momenta defined in Eq. JBJ) for the eikonal velocities of the quark 
pair. 

The long-distance evolution of the pair from octet to singlet color configurations is given by 
the infrared factor E^ 1 ^ of Eq. (jHj), which can be represented as a matrix element. This matrix 
element is given in the notation of Eq. (jSJ) by 

£ {8 ^ 1) (P l ,P 2 ,s) = £>| W^(0)] tt [T^] jk [*gt (0 )l $^t ( o) e6 \N) 



N 



IJ J JJi L \KI 

(N\M A \0) ba *g(0) r \T^] MN *g(0) ^ |0>. (12) 



LM 



NL 



Here we have exhibited all color indices: those in adjoint representation by a, b ... , and those 
in the fundamental representation by I, J . . . , to indicate the trace structure, which imposes a 
color singlet configuration on the quark pair in the final state. In Eq. (112)) and below, overall 
time-ordering of the field operators is understood in the amplitude, and anti-time ordering in its 
complex conjugate. 

The operators $^ <? - ) and are the ordered exponentials that represent the quark and anti- 
quark. The quark (antiquark) ordered exponential has the same (opposite) ordering compared 
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to time ordering. To be specific, we represent normal (reverse) matrix ordering by V (V), and 
define 



$g(0) = Pexp 
$g(0) = Vexp 



coo 

-w 



'0 

poo 

w / 

JO 



poo 

/ dXP 1 -A ( - q \P 1 X) 
Jo 

dXP 2 - A {q \P 2 X) 



(13) 



For classical fields, $®(0) is the hermitian conjugate of $n\o). The matrix A„ = ^2 a Ta 9 ^ A v ^ a 
is the gauge field operator in the quark fundamental representation. 

The matrix element in Eq. ()12j) is equal to unity when q = 0. In Ref. we expanded Eq. 
(H2J) to second order in q, and found an expression in terms of field strength operators, 

poo 

et*% + Q,P-9,e) = E/ d\'\' (0\<S>\ A) \0) M ,<S>l A \\%AP^ u ^Am\N) 

M JO 

poo 

(N\ / dX X ^ A \0) bd WF V(t>a {Xp) ] & p A \X) ad |0) . (14) 
Jo 



N 

X 



In this expression, the momenta and q^ are taken to be dimensionless, scaled by the heavy 
quark mass m. We notice that in the heavy quark rest frame, p = (l/2m)(Pi + -P2)rest = ^o, the 
relevant operator is precisely the chromo-electric field F^q, and the matrix elements describe an 
electric dipole transition. 

The basic result of Ref. f)] was to identify an infrared pole in £^~* l \ associated with the 
exchange of gluons between the heavy quark pair and the eikonal source The presence of such 
a pole showed first that the infrared behavior of the fragmentation function is not summarized by 
"topologically factorized" diagrams alone in an arbitrary gauge ■ The specific form, however, 
depends on the relative velocity of the pair only. In terms of the center-of-mass velocity of the 
heavy quark and antiquark, defined as in Eq. © above, the result is 1 

4^\v) = ^{Nt-l)al±-v\ (15) 

The prefactor here is the same color factor given in Eq. Q above. The crucial point is that 
the pole due to topologically non-factored diagrams is independent of the direction of the light- 
like vector / M . In Ref. jH] we showed that this implies that the gauge-completion of the NRQCD 
matrix element is adequate to match this infrared structure in production processes with arbitrary 
numbers of recoiling jets. This verified the NRQCD factorization of Eq. (j2J) at NNLO in soft 
gluon corrections for octet operators. We will show here that this result generalizes to all orders 
in v for the full infrared factor S^ 1 ^. 

In our calculation below, we will keep q, or equivalently the relative velocity v , finite and 
nonzero, and evaluate Eq. (JT2jl directly to NNLO (a,), without an expansion in velocity. Re- 
stricting ourselves to this order in a s , an expansion in v will be describable at any order in terms 

1 The coefficient of v 2 here is four times larger than in Ref. [0], where the velocity was defined by 2q/m, as 
noted in connection with Eq. © above. 
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(a) 




(b) 

Figure 1: Examples of NNLO diagrams for the quarkonium production, a) Three-gluon diagram 
b) QED-like diagrams. The projection of the quark pair onto a color singlet in the final state is 
understood. 

of derivatives of the electric field strength. In this sense our calculation identifies the infrared 
pole in the sum of electric multipole transitions generated by NNLO. Because of the eikonal 
approximation, no information on spin-dependence is included in the calculation, and the con- 
sequences for matching are similar to the v 2 case. We will once again find an infrared pole to 
any order in v 2 , and we will once again find that the pole is independent of the direction of 
the vector Z M . This surprising result shows that matching to gauge-completed NRQCD matrix 
elements is not limited to lowest order in v 2 , but to this order in a s and for this class of electric 
dipole transitions, is true to all orders in the relative velocity v. 

3 Finite-?; diagrams 

3.1 Ladder and three-gluon diagrams 

Our goal is to calculate the non-canceling infrared pole term in Eq. ([12)1. the eikonal infrared 
factor at finite relative velocity v. As in Ref. 0, we need only consider diagrams that are not 
topologically factorized. As indicated above, we take the quark momentum as Pi, the antiquark 
as P2, and the momentum of the gauge line as I, with I 2 = 0. The gauge line may be thought of 
as part of a fragmentation function, or as representing an approximation to a recoiling jet in a 
hard-scattering cross section. 

The relevant diagrams are illustrated in Fig. ^ with quark-antiquark eikonal lines on top, 
and the gluon eikonal below. In each case, the vertical line indicates the color-singlet heavy 
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quark pair, and the sum over other possible final states is understood. The left- and right-most 
vertices in each diagram couple the light-like color octet eikonal on the bottom of the diagram 
with the color octet projection of the eikonals on the top that represent the heavy quark pair, in 
the fundamental representations. 

The quark lines radiate soft gluons to transform themselves from color octet to color singlet, 
as in the classic description of the color octet mechanism. At NNLO, however, these soft gluons 
can scatter from the eikonal gluon in the I direction. To this order, the diagrams either contain 
the three-gluon coupling as in Fig. or are QED-like as in Fig. QJd. Of these, only the former 
give rise to noncanceling infrared poles, generalizing the result of Ref. [Oj. 2 The color factor in 
Eq. (|TU|) is thus found entirely from Fig. Ek- This factor is the product of two simple traces in 
fundamental representation, and the color tensors from two three-gluon vertices, 



The two traces reflect the projection of the quark- ant iquark eikonal pair onto the singlet in 
the final state that is built into the matrix element of Eq. (jl2j) . In the traces, one generator 
is associated with the operator in the corresponding matrix element and the other with the 
corresponding quark-gluon vertex shown in Fig. We now turn to the calculation of the 
velocity dependence of E^ 1 ^. 

3.2 Velocity-dependence of £'( 8 ^ 1 ) 

We will give a detailed discussion of the cut diagrams illustrated in Fig. EJ In these four cut 
diagrams, a soft gluon is emitted by one member of the (eikonal) heavy quark pair in the am- 
plitude and is absorbed by one in the complex conjugate amplitude. This soft gluon rescatters 
from the spectator eikonal gluon of momentum / in either the amplitude or complex conjugate. 
For each of these diagrams the leading order v 2 expansion of this contribution to £2, Eq. (fT^j) . 
was calculated in detail in Ref. || (where it was referred to as diagram IIIA). Here we extend 
this calculation to all orders in v by directly calculating the full velocity dependence of E^ 8 ^ 1 ^. 

All of the diagrams in Fig. 2 have the same color factor (the same as in Eq. (j!6|) above) when 
the quark pair is projected onto a color singlet in the final state. We will suppress the common 
color factor in our calculation below, and denote the sum of the relevant NNLO diagrams as 
j(8^i) _ £(8->i) 1 p co[or . We label the diagrams of Fig. |2J whose sum gives T$~* l \ by 



where the first superscript, Pj, on the right identifies the gluon coupling in the amplitude, to the 
quark, P 1; or antiquark, P 2 , and the second, Pj, the gluon coupling in the complex conjugate. 

2 It is worth noting that the color factor of the left-most diagram in Fig. vanishes in the matrix element of 
Eq. CTJ. 



color 



Tr [Tj*> ] Tr f eai f ihg , 




(16) 




(17) 
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As in Ref. j^j, we fix the momentum of the eikonal I in the minus light cone direction, 



F = 6 



p-- ■ 



The momenta P\ and P 2 remain arbitrary. 







Figure 2: a) The four relevant cut diagrams with three-gluon vertices. The Pi, P 2 and / lines are 
eikonal. The full contribution of each diagram to X^ 8_>1 ^ is found by summing over other cuts that 
include the singlet pair, b) An example of a two-gluon final state that cancels the (k 2 — ki) 2 = 
pole, as discussed in the text. As in Fig. 1, the projection of the quark pair onto a color singlet 
in the final state is understood. 



3.2.1 The P\P\ gluon rescattering diagram 

We begin with the diagram X PlPl , Fig. EK, which describes the interference between the order g\ 
rescattering of a gluon that is emitted by the heavy quark (Pi) eikonal line and the lowest-order 
process in which it is emitted by the same line. 

Written out with the momentum structure of the eikonal vertices in Eq. (|11J1 shown explicitly, 
this diagram is given by 

I P ^(P,q,l) = f 2n5(k 2 2 ) [{h _ k * )2 - - } j^-j 

x (-g/jfV^xlki, k 2 - k h -k 2 ] ) 



x 



A 



-i(-gfi £ 



[Pi-h+ie] [Pi-h-ie] [I ■ (k 2 - h) + ie] 



(19) 
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where the three-gluon vertex (with all momenta flowing in) is 



^Ml,/*2,A»3 [<7l> 111 Is] — (<?1 — ^fls 9nxll2 + (?2 — ^J/i! 9n2fl3 + (?3 — Ql)(l2 9/J.3V1 



(20) 



As in Eq. (jHJ), P and g are the total and relative momenta of the pair. As in J^ 1 ), Eq. (TJJ), Pi 
and P2 are the momenta of the heavy quark and antiquark, respectively, and I is the momentum 
of the lightlike eikonal. As indicated in the diagram, we take k\ as the momentum of the soft 
gluon emitted in the amplitude, to the left of the cut in Fig. and k 2 as the momentum of the 
soft gluon flowing to the right in the figure. 3 

Following the procedure of Ref . , we find it useful to integrate first the minus components of 
k\ and k 2 . Performing the k^ integration by using the mass-shell delta function, and evaluating 
the numerator factors, we find 



4 4e 

9 K 



d D - 2 k\ 



dkf 



Pi 



k T 



2 /lo 



(27i) 2D ~ 1 

+ fcr) + P{{kt + kt) - pi ■ (# + k{ 
1 




kj + ie] 



p i k 2 2 I p-u+ _ pT . L.T 

2k+ + 1 2 1 2 

1 



IE 



T 2 



2(ki ^J)(^i 7^f+) (kj k 2 ) 



ie 



[2k+K ~ kf + ie 



h- + ELK 



^1 +ie 



(21) 



Below we will generally suppress limits on the loop (hi) and real-gluon momentum (k 2 ) integrals, 
except where they are necessary for the argument. The k\ integrals are unbounded, while the 
k 2 integrals are understood to be cut off at the order of the quark mass. These upper limits are 
set by the phase space cutoff, denoted \i here, and discussed above in connection with the basic 
matrix element of Eq. (JSJ • 

The next step is the loop integral k^ in Eq. (|21|). which we perform by contour integration. 
The integrand of (J21)) has three poles in fcj", 



K [kf] 

k l '[(fci-fc) 2 ] 
K[Pi-ki] 



T2 



ie 



2k+ 



k 



T2 



2 , (^l ^ 2 ) 2 l f 



+ 



2/^2 2^/c-^ ^2 



IE 



(22) 



labeled according to the denominator that vanishes. 

3 In the discussion that follows, we will compute diagrams with loops in the amplitude, in contrast to the 
complex conjugate amplitude as in Ref. Since our result is real, the analysis is otherwise completely equivalent. 
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The pattern of poles encountered here is similar to those for the order-f 2 calculation of Ref. 
jU], and the role of each term is similar. When we close the fcjf contour in the lower half-plane, 
we pick up the k^^ pole when ki > 0, the fcfKfci-fe) 2 ] P°l e when kf > ki and the ki^.^] pole 
for all values of kf . The contribution from the k^^ pole vanishes, because the resulting integral 
is antisymmetric under the exchange of the remaining components of ki and k 2 - The ^flOi-^) 2 ] 
contribution, on the other hand, cancels against the corresponding cut in which the gluon with 
momentum k\ — k 2 appears in the final state, as in Fig. 2b. 

This leaves us with the third, fcf [p^ki], pole only, whose calculation we describe in detail below 
for Z PlPl (P, q, I) and J PlP2 (P, q, I). We note that in the sum of J PlP2 (P, q, I) and J P2Pl (P, q, I), 
the symmetry argument for the cancellation of the k^^ pole continues to apply. 

In summary, to derive the infrared contribution to X PlPl (P, q, I) from Fig. we close the k^ 
integration in Eq. (|2*T|) in the lower half-plane at the pole fc^^.^j, and find 

J**(P,g,i) = (-27ri)(i)J^ J d D ~*ti[ J dkt J d D -V 2 J 
'k T2 



dk+ 
2ki 



p;- ( ^ ( ^# - ^f] 1 + PiiK + kt) - p p . {k T 2 + kD - ^{kt + kt) 



2k 2 
1 



pi- 



pi 



i 



[ki -k+ - ie\ 



+ P,-ki - Pf ■ Af - ie 



x 



x 



2W-* 2 + )((^-^)-$) -W-*D 2 + 



1 



(23) 



In this expression we have two (P-2 dimensional) transverse and two plus integrals remaining. 
We begin by applying a Feynman parametrization to the final two denominators, 



l P ^(P,q,l) 



X 



X 



Pt 



4 4e 

9 1* 



(2vr 



,2D-2 



2 /uo 



+ 



1 



rfx / d D ~ 2 fc P / dfc? 



[ki-kt-ie] 
1 

[fc^ 2 + UPi) - it} 2 



dki J 


dk+ 
2ki 


P{ki 


))* 


Pt 






+ ki 



^f + P 1 -ki-P[.k 2 T - i e 



d D - 2 k[ T 
Pi (ki + ki 



(24) 
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where to complete the square we have shifted to k[ = kf — K (Pi), with 



K 1 (Pi) =x[ki + 



(kf 



Pf 



(25) 



In adddition, in the final denominator of (|24j) we introduce the function 

2 



L(Pi) = x 



L^2 



X 



k^ 

l.T K 2 D T 
2 p+ 



Pi 



+ 



^2 ^1 p2 

xPf 1 



(26) 



Here and below, to simplify the notation we suppress the k\ and /c2 dependence in L{Pf) and 
A' '(/':)• 

We now perform the integration of Eq. (}2"4"|) in D=4-2e dimensions, obtaining 



71 



X 



l-£ 



r(l + e) 



4 4s 



(2< 



2D-2 



jf ^ y dkt j d°~ 2 ki j 



dkt 
2 A?q 



A ^+ + Pi ^2 — A ' 



P? 



(4 + kf) 



X- 



[^1 ^2 
1 



+ Pftf - PI ■ kl - ie 



x 



[L(Pi 



16 



1+e 



(27) 



Simplifying this expression algebraically, we put it into a form that will facilitate the combination 
of diagrams below, 



X^(P,q,l) 



71 



1-e 



T(l + e) 



4 4e 

9 f* 



(27T] 



2D-2 



dx Id 



]D-2 iT 



X 



[kf - kf - ie] [LiP^-ie] 1 - 



1 - 



dki 



p+ (kf + kf 



dkf 



P+k r2 



Pi ^2 -^1 ' ^2 ^ C 



• (28) 



Of particular interest in this expression is the first term, 1, in square brackets. This term has up 
to four poles in dimensional regularization, corresponding to momentum configurations in which 
both momenta k\ and &2 vanish and are collinear to the light-like eikonal line I in Fig. El We will 
see that, as in Ref. 0, these singularities cancel, leaving only a single real 1/e (and imaginary 
1/e 2 ) pole in dimensional regularization. 



3.2.2 The P1P2 diagram and collinear cancellation 

Before continuing with the integrals, we turn our attention to the third diagram in Fig. |2t, which 
describes interference between gluon (ki) rescatterring after emission from the heavy quark (Pf) 
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line in the amplitude with gluon (k 2 ) emission from the antiquark (P 2 ) line in the complex 
conjugate amplitude. We denote this diagram by Z PlP2 (P, q, I). 

The momentum-space integral for the P1P2 rescattering diagram, I PlP2 (P,q,l) is 

xP^p^v^kxM-K-h] 

1 

X [Pi • h + ie] [P 2 ■ k 2 - it] [I ■ (k 2 - h) + ie] [{h - k 2 y + it] [k\ + ie\ ' 

(29) 

This diagram has an overall (—1) relative to I PlPl (P, q, I), associated with the connection of one 
gluon to the P 2 eikonal line, which is in the antiquark representation. 

Performing on Z PlP2 (P, q, I) the same steps as above for the k~ and k\ integrals in Z PlPl (P, q, I), 
we can put this integral into a form analogous to Eq. f)28j) . although slightly more complex, 

X p ^{P,q,l) = n^T(l + e) j\x j j g| j dkt 

x ^ 1 [/_ 

[k+ -kt-te] [L(Pi) -ie} 1+£ [ 



^(kt + H) + T?(kf- k+KP 2 + Pr - P?P 2 ~) + 2 P 2 + 


> 2 T 


■ [^) - kl] 




p 2 k 2 2 I p-U+ pT ]„T 
2k + + ^2 ft 2 -*2 fi '2 ' e 





(30) 

Combining the expressions for Z PlPl (P, g,Z) and I PlP2 (P,q,l), Eqs. (J28)) and (|3U|) . we see that, 
as anticipated above, the collinear-singular terms (the Is) in the square brackets cancel. 

Our goal now is to evaluate the infrared poles from the expressions in Eqs. (|2~%j) and (JHOJ). 
For this purpose we need to perform the kf integration, which ranges from -00 to +00. Noting 
that the numerators are at most linear in kf, we rewrite factors of kf in the numerator as 
(kf — k£) + k£. The (kf — k£) term then cancels the corresponding denominator in Eq. (J3U|) . 
We next reorganize the contributions of the Pi Pi and P1P2 diagrams to Z( 8_>1 ) as 

jAft + jftft = (/ift + jPxP.) + ( K PyPi + JCP^) , (31) 

where, for example, J PlP2 is obtained from I PlP2 by canceling the (fc^ — k% ) factors in the 
numerator and denominator, and )C PlP2 represents the remaining terms in I PlP2 . We now turn 
to the identification of infrared poles in these expressions, using slightly different procedures in 
the two cases. 
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3.3 The IR pole from the J terms 

We begin with the terms that lack the pole in kf — k 2 . From Eq. (J30|) we find 



J p ^(P,q,l) + J p ^(P,q,l) 



71 



1-e 



T(1+e) 



4 4e 

g v 



(2vr; 



2D-2 



dx I dk+ I d D ~ 2 kl 



X 



dk+ 



P1P2 I 2 ( p+ p- p+ p- 
P+P+ "+" p+p+ \ r 2 r i ~ r l r 2 



P T 



IL(Pi)-^} 1 - 
-(P 2 - A) • 



p+ z/i 2 



P? 



^2 



P~T 



IE 



(32) 



In these terms, the remaining (quadratic) /c^-dependence is in the L(P\) denominator. The 
integral is elementary, and we find 



J p ^(P,q,l) + J p ^(P,q,l) = -2 1+2£ {-l-ie)- 1 ' 2 ^ 5(1/2, 1/2 + e) 7T A -T(1 



X 



4 4s 

07* 



p2 



da; 



™l+2e I i +l+2e 

o x Jo kZ 



l-ej 



+ + 



p+ p+ 1 p+ p- 



■(P+P,- -P+P 2 ) + 2 



p? 



EL 



Pi+ ) + Pi+2 



-1 l+2e 



P 2 \2 1 P 2 2 
^2 +J + ^ 



Pi) 



2 



(33) 



The leading behavior of this expression is a purely imaginary double pole in e, from the lower 
limits of the x and k 2 integrals. The remaining transverse integration is both infrared and 
ultraviolet finite, and real, for e — > 0, so that an expansion in e of this integral can give rise to 
only imaginary single poles, which vanish in the cross section. The overall factor of (— 1 — ie)~ e , 
however, can convert an imaginary double pole to a real, single infrared pole. To isolate the 
residue of this pole, we need only evaluate the k\ transverse integral at e — 0. An important 
point is that, because we are interested only in noncancelling infrared poles, we may extend the 
upper limit of the k 2 integral to infinity. We can do this because at finite k 2 > fi, only the k\ 
line can produce an infrared pole, at one loop. But, as discussed in jS], for example, all one-loop 
infrared divergences factorize in the sense of NRQCD. 

A simple change of variables, yx = k 2 /k 2 , simplifies the transverse integration, evaluated at 
e = as just described, 



d 2 k T 2 
{ktf 



1 



fe 2 



pT 



+ 



p? 



Pi 



^2 




^2 



d 2 y T 



p2 



[y 2 T 



i] 



{y T + a T ) 2 + ^ i 



(34) 
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where we define 



a? 



pT pT 
r l r 2 



'P+ 



2 ■ 



Pi P2J V ^1 

Next, introducing another Feynman parametrization, w for and integrating over y^ we get 

1 



(35) 



d y T 



[y 2 T + 1] 



7r / dw- 



p? p? 



1 — w — w 2 aip + w(a^ + 7^- — 

^2 



(36) 



After integrating over w we then find, for the original k\ integral, 

d 2 k\ 1 n 



(kty 



k+ p+ J p+< 



UT pT\ 2 p2 

" ■ > •* 9 l j 2 



In 



a + c 



a — c 



(37) 



where 



2Pi ■ P 2 

P+P+ 

r \ r 2 



a \ 1 



p2 p2 
-'l ^2 

(Pi ■ P2 



\2 ■ 



(3? 



We are now ready to isolate the real single pole, as indicated above. The poles of the remaining 
x and k£ integrations in Eq. ([315)) are found by using 1/ f 1+a = —6(f) /a, and for the e-dependent 
phase, the expansion 

1 ± me + 0(e 2 ) + ... 



(-l±ie 

gives the corresponding £ term. We then find from Eq 
(which are determined by simple substitution), 

2Re[j p ^ + J**+J** + J p * p >] 

1 + 



(39) 



adding the two remaining diagrams 



Te 



2 r 



p+p+ 

r l r 2 



Pi ■ p 2 / 2/1 V 1 - /i 



(40) 



where AP 2 and h depend on the momenta of the pair as 



AP' 



h 



A _ Ei 

P + p+ 

r l r 2 

p2 p2 
f_l r 2 

^ (Pl-P2f 



(41) 
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It is easy to see that h is proportional to the velocity of the heavy quarks in the center of mass 
system. We will give an explicit expression for its full f-dependence below, after identifying the 
poles of the /C terms. Taken in isolation, however, the sum of the J terms is manifestly not 
independent of the choice of Z M , because of the explicit dependence on the factors P^ = Pi ■ i in 
fJ4f)|) . 4 Independence from the direction of the light-like Wilson line used to define the underlying 
matrix element will, however, emerge below, once we find the infrared poles of the K, terms. 



3.4 The IR pole from the /C terms 

Now we evaluate the remaining terms in Eq. (jSHJ), for which the factor k± — k 2 remains in the 
denominator, such as K PlP ' 2 = 1 PlP2 — J PlP2 . In these cases, the virtual kf dependence in the 
denominators is no longer quadratic, and the integral requires a slightly more elaborate approach. 
Specifically, the contributions from the P1P2 and P\P\ diagrams are given by 



K P ^(P,q,l)+}C P ^(P, q ,l) 



x 



dx 



dkt 



4 4e 

9 V 

(2tt) 2D - 2 



2 /un 



dkt I d D ~ 2 k T 



[kf -ft- ie] [L(Pi) - ie] 



l+e 



^(2kt)+2P+ 






■ [G T (Pi) - kl] 




' p+hT 2 

L 2 4 + 2 


^2 


-PI 


■ k\ — ie 





where 



(P2 - Pi) ■ 



G T (P 1 )=xkJ+(l-x)kt%. 



(42) 



(43) 



At this stage we employ another Feynman parametrization (y) to organize and perform the k 2 
integration. As for the J terms, we may extend the upper limit of the transverse integral to 
infinity, with the result 



c/ 4 (47r/i 2 ) 2£ r(l + 2e) / dx j dyy e / dk\ 



32tH 



dki 



x 



Pi-P 2 

P+P+ 



+ AP 2 (1 -x)(l -y) 



k 2 [k l k 2 ie] 
1 



(M-ie) 1+£ (N + ie) 1+2£ 



(P2 - Pi) 



(44) 



4 It is interesting to observe that the overall result in (|40|l is invariant under independent rescalings of Pi, P2 
and I, as expected for the eikonal approximation. 
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where we define 



x 



N = kn 



-y(l-y)AP 2 + 



V / k+P? 



-{x- x/y) 



Pi 



(45) 



with AP 2 defined in Eq. (gTJ). In Eqs. (02) and (jUJ), we have exhibited /C PlP2 , from the third 
figure in the diagram Fig 2a. The other term may be found as indicated, simply by replacing P 2 
everywhere by Pi, which leads to a number of simplifications. For example, AP 2 is replaced by 
zero in these terms. 

Equation (03} can be reorganized slightly to yield 



32tt 4 



x 



/(47r/i i ) 2£ r(l + 2e) / dx / dyy 



1 



o 



-l+4e 



(x(z — x) — ie) 1+£ (y/x) l+2e 



dz 



z — 1 — ie 



x 



1 - y)AP^ + x(l/y -l)^% + z + 

9 1 



1 l+2e 



X 



f_l 

A 2 



- (P 2 - Pi) 



Pi 

Pi ■ Pi 

P+P+ 

r l r 2 



AP\\-x){\-y) 



16 



(46) 



where we define z = /kj- The variable z appears in three denominators, one of which is a 
simple pole at z — 1 in the upper half plane. To perform the z integral in Eq. (|3EJ), we combine 
the remaining two denominators by introducing another Feynman parameter, y', 



[x[z — x) — ie) 



-x(l - y)AP 2 4r + x(l/y - l)£pr^f + z + ie 



p+ z Pi 



l+2e 



T(2 + 3e) 



x 1+e r(l + e)r(l + 2e) 
dy' 



y' s (l-y 



>\2e, 



-1+ie) 



-l-2e 



(2y> - l)z + x(l - y>) ((1 - y)AP^ 



1 2+3e 



(47) 



Notice that before the combination of these two ^-dependent denominators with fractional powers, 
we factor —1+ie from the second of the two, so that the ie on the right-hand side of (07} has 
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a definite sign: —it. It is perhaps worth noting further how we determine this prescription. As 
z varies from -oo to +00, the second (and first) denominator in ()47|) can vanish, but for any 
variable denominator, denoted W, we have 

W + ie = W if W>0 

W + ie= (-l + ie)(-W -ie) if W < 0. (48) 

The factor (— 1 + ie) £ = 1 + ine + . . . in this expression will be particularly important in the 
calculation below, because as above it will convert an imaginary double pole in e to a real single 
pole. 

Although Eq. (}4"Tj) is a bit complicated overall, the combined denominator on the right-hand 
side gives a branch cut in the z-plane that is in the lower (upper) half-plane for < y' < 1/2 
(1/2 < y' < 1). Inserting the expression in Eq. (|4"TJ) back into Eq. (jUj) , we then easily perform 
the z integral. For y' > 1/2 we can complete the z contour in the lower half plane without 
enclosing any singularities, and the z integral vanishes. For y' < 1/2 we close in the upper half 
plane and pick up the simple pole at z = 1. 

The result of this procedure is 

- *>--* r (2 + as) MA _,^ f dki r ixx , r dyyl+ * 



16vr 3 T(l + a 

1/2 

dy'y'%l-y') 2£ x 







- 9\^ 2 r j 


f dk+ 


z,+ 1+4e 


/ _l2\ 1+2e 

(?) 


Pi ■ P2 




p+p+ 


1 




-3/0/(^1,^2 


) + xyy' 







+ AP 2 (1 -x)(l -y) 



(1 - 2y')y + x(l -y)(l- y')f(Pi, P 2 ) + xyy' + te ] 2+3e 
{Pi -> A) • (49) 



where 



p+ 2 p 2 p+ 2 
/(^^) = -^-^-Z/AP 2 -^. (50) 
P 2 n 1 

In the above equation, (— 1 — ze)^ 1 ^ 2 is obtained by multiplying the overall factor (—1 + ze) _1_2£ , 
from Eq. ([47)1 . together with (— 1 — ze)~ 2_3e , from the denominator on the right-hand side of (jUj) , 
after the z integration. Although the result has the same number of integrations, the infrared 
behavior of the y' integral is straightforward to analyze, which will simplify our determination of 
the infrared poles. 

We can check the infrared behavior of the y' integration by expanding in e, 

y"(l- y'f- = 1 + £ [\ny' + 2 ln(l - y')} + .... (51) 
The term in square brackets is finite at the endpoint y' = 1/2. Because the denominator is 

1 /2 

finite at y' = 0, the behavior at y' — is then no worse than logarithmic, J Q dy'lny', and the 
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y' integral is finite and real for e = 0. Thus, to identify the real infrared pole, we can replace 
. — y') 2e by unity and integrate over y' in (}4*9"j) to obtain an expression that includes the full 



real single pole, 

= i (—1 — ie 



-1-e 



T(2 + 3e) 



l + 3e) r(l + e) 16tt 



^ 4 (4tt/x 



2\2e 



dki 



Pt 



Pi -Pi 



r l V r l r 2 



+ AP 2 (1 -x)(l -y) 



x 



_l+4e 



1 



dyy 



l+2e 



[-2y-a;(l- ? /)/(P 1 ,P 2 )+zy] 
1 



[§z(l - y)f(P 1 , P 2 ) + \xyY+* [y + x(l ~ y)f(Pi, P 2 )] 1+3£ J 



- (Pz —> Pi] 



(52) 



Here we note that terms of order e from the y integration can also produce at worst imaginary 
single poles in e. 

The steps that isolate the real infrared pole for the /C's follow the procedure of the previous 
subsection. The x and y integrals in Eq. (|52j) integrals are real, and we must identify a 1/e 
pole from these integrals. The resulting double pole then combines with the me term from the 
expansion of the prefactor (—1 — ie)~ 1 ~ e to give a real pole term that survives in the full sum of 
diagrams 

It is easy to check that the x and y integrals of the second fraction in square brackets in 
Eq. ()52j) are both real and finite for e — > 0. We may thus limit our analysis to the first fraction, 
which is proportional to l/x 1+3e . The x — >• limit thus generates another infrared pole, which 
can be isolated by using the distribution identity, l/x 1+a = —5(x)/a, while the y integral can be 
performed at e = 0. 

For the remaining analysis, it is convenient to combine Eq. (|52|) with the corresponding results 
from the other two diagrams in Fig. |2K, which can be found by substitution. After performing 
the x integration in (j52j) we obtain, suppressing the finite part, 



ii-l-ie)- 1 - 6 



r(2 + 3e) 



2 3e g 4 (4tt/z 



32tt 3 £ 

2 \ l+2e 



Pi 



PI 



(l + 3e) r(l + e 
Pi- Pi 



2\2e 



dk+ 

+ l+4e 



dyy 



2s. 



+ (Pi 



I P 1 P 2 
~ (P2- 



AP 2 (1 - y) 



Pi) 



(Pi 



J 



[l-y)f(P 1 ,P 2 )] 



l+3e 



(53) 



Taking into account the k£ integral, we see that this expression has an overall imaginary double 
pole which, however, cancels when the complex conjugate diagrams are combined. 

The leading, real pole in e can be obtained from the above equation by setting e — > in 
the y integral, which then reduces to the sum of two elementary integrations. The first has a 
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y-independent numerator in the integrand 
dy 



i pfjpl 



[y+(l-y)f(P 1 ,P 2 )} 



dy 



[y+(l-y)f(P 2 ,P 1 )} c 



a + c 



a — c 



(54) 



with a and c defined in Eq. 
given by 



PfjPl 



o dyy h + (i-v)f{Pi,p*)] 



The additional term is linear in y in the numerator, and is 
d(Px,P 2 ) 



2c 



In 



a + c 



a — c 



2AP 2 



In 



p2 p+2 

p2 p+ 2 
-"2 -"1 



(55) 



where a and c are defined in (J38|) . AP 2 in (|41j) and 

1 P 2 



d(P!,P 2 ) = 1 



P 2 



AP 2V P+ 2 P 2 



.21- 



(56) 



Although the integral of Eq. ()55j) is not symmetric in P\ and P 2 , once we add the four diagrams 
we obtain symmetric result, equal to (J51j) . 



dy y 



PfjPl 



+ 



y + (1 - y)/(P!, P 2 ) j/ + (1 - y)/(P 2 , ^ 



iln 

c 



a + c 



a — c 



Combining the results above in Eq. (jEBJ), and using 



■l±ie) 



we obtain for the infrared pole, 



,,±tTTE 



1 ± ins + C(£ 2 ) + . . . , 



(57) 



(5? 



2 Re f" K PlPa + /C PaPl + /C PlPl + /C P2P2 1 = ^ 
L J Ae 



4e 



2 + 



p+p+ 

*i ^ 2 -AP 2 I In 



P-P 2 



1 

2ft 



c \ a — c 
1 + ft 



1 - ft 



(59) 



with ft defined in Eq. (141 jl . Again, we find a dependence on the choice of vector I, which, however, 
cancels in the sum of the J7"s and )Cs. 

3.5 The IR pole in I 8 ^ 1 

Now adding Eqs. flU and $2$ we find from Eq. (HU1) 



T^ 1 ) = 2 Re [ J PlPa + J PaPl + J PlPl + J P2P2 ] 



As 



2ft 



In 



1 + ft 
1 - ft 



(60) 
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with the function h defined in Eq. (|4ip. This is the general expression for the pole part of Z^ 1 ). 
A short calculation shows that the lowest-order expansion of Eq. (160)1 reproduces the result of Ref. 
[Oj, Eq. |T5|) given above for the infrared pole at order v 2 (Electric dipole). As anticipated above, 
our result is independent of the direction of the light-like vector V 1 that defines the integrals. 

To derive an equivalent form in terms of the relative velocity, we use P = Pi + P2 and 
1q = P\ — P 2 and recall that we identify the relative velocity with the heavy quark velocity in 
the pair center of mass: v = q/E*, where 2E* is the total energy of the heavy quark pair in this 
frame. In these terms, we may replace the quantity h of Eq. (jSOjl with an explicit function of 
velocity, 



2<8-i) = ^± 
Ae 



1 + 

2/(N) ~LWOTJ 



In 



(61) 



where v = \v\, and the function f(v) is given by 



/(*) 



2v 



1 + v 2 



(62) 



Equation (|6Tjl is the general result at NNLO for the infrared pole. As Eqs. (joTj) and (jS2]) show, 
the infrared term is independent of the eikonal momentum I for finite v. Hence the infrared pole 
structure at NNLO in a s is consistent with factorization, which at this order in a s is valid to all 
orders in the relative velocity v in the heavy quarkonium system. 



4 Conclusions 

Equations (|61[) and (|62|) provides a remarkably compact expression for the single pole in the 
infrared factor Z( 8_>1 ) of Eq. (fT2*J) at NNLO for finite relative velocity v, or equivalently, expanded 
to all orders in v. This result is independent of the direction Z M of the octet Wilson line, and hence 
is consistent with NRQCD factorization, as discussed in connection with Eq. (|Sjl. As shown in 
Ref. 0, the absence of /-dependence in gauge-completed matrix elements enables them to match 
all NNLO infrared divergent corrections to multijet cross sections. All the arguments in that 
reference for matching at order v 2 apply here to all orders in v. 

Although limited to NNLO, our result suggests that the decoupling of light parton dynamics 
from heavy quark pair production is robust in perturbation theory at the level of infrared diver- 
gences. Evidently, at NNLO, while a light-like, energetic parton can resolve the color structure 
of a heavy quark pair, it does so in a way that is independent of the direction of the relative 
motion of the pair. This suggests to us that the results derived here may generalize to higher 
orders in soft gluon exchange. 

In closing we note that there are, of course, many terms in a general NRQCD velocity ex- 
pansion that are not given by the eikonal approximation, in particular those that deal with spin. 
The eikonal approximation, however, does reproduce infrared divergences as they appear in per- 
turbative calculations of heavy quark production. Clearly, an extension of this analysis to higher 
orders in the coupling will be of interest. 
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